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Abstract Energetic gamma rays (GeV to TeV photon energy) have been detected to- 
ward several supernova remnants (SNR) that are associated with molecular clouds. 
If the gamma rays are produced mainly by hadronic processes rather than lep- 
tonic processes like bremsstrahlung, then the flux of energetic cosmic ray nuclei 
(> 1 GeV) required to produce the gamma rays can be inferred at the site where 
the particles are accelerated in SNR shocks. It is of great interest to understand the 
acceleration of the cosmic rays of lower energy (< 1 GeV) that accompany the en- 
ergetic component. These particles of lower energy are most effective in ionizing 
interstellar gas, which leaves an observable imprint on the interstellar ion chemistry. 
A correlation of energetic gamma radiation with enhanced interstellar ionization can 
thus be used to support the hadronic origin of the gamma rays and to constrain the 
acceleration of ionizing cosmic rays in SNR. Using observational gamma ray data, 
the primary cosmic ray proton spectrum can be modeled for E > 1 GeV, and care- 
ful extrapolation of the spectrum to lower energies offers a method to calculate the 
ionization rate of the molecular cloud. 
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1 Introduction 

The origin of cosmic rays (CRs) is one of the central questions in modern astro- 
physics. In particular, for gamma rays the formation processes are still a matter 
of debate due to ambiguity of the observations. The detection of the emission of 
gamma rays in GeV and TeV towards molecular clouds in the direct vicinity of su- 
pernova remnants (e.g. HI) has further fueled the question about the origin of this 
radiation. In principle, there are three different scenarios by which these particular 
gamma rays can be formed: bremsstrahlung or inverse Compton scattering of elec- 
trons in a leptonic scenario or the formation and - in turn - decay of neutral pions 
from scattering of CR protons on ambient protons in a hadronic scenario. In most 
cases, it is not possible to determine which process actually is at work based on the 
gamma ray spectrum alone. Therefore, correlation studies need to be performed. In 
the hadronic scenario, CR protons of energies E > 1 GeV form neutral pions by 
deep inelastic scattering on ambient matter, which in turn decay into two gamma 
rays. It is to be expected that there are also lower energy protons, which are most 
effective in ionizing hydrogen. Using molecular clouds as tracers for these ionizing 
lower energy CR protons in spatial correlation with GeV to TeV gamma radiation 
could offer a strong hint at hadronic origin of the gamma rays, since CR electrons do 
not penetrate molecular clouds as deeply as CR protons, resulting in a different ion- 
ization profile. The ionized molecular hydrogen, Hj, triggers a chemical network, 
forming about 15 to 20 molecular ions, which are formed in ro-vibrationally excited 
states (see e.g. JIBU TI). If there is a sufficient amount of these ions formed, i.e. the 
ionization rate is high enough, the relaxation of these ions results in characteristic 
line emission, which can be detected as tracer for enhanced ionization [6 1. The basic 
idea is sketched in figure [T] 



2 Deriving the CR proton spectrum 

To calculate the ionization rate of the molecular cloud, the primary CR proton spec- 
trum at the source is required. To obtain this, observational data about the gamma 
radiation from instruments as e.g. FermiLAT can be used. To generate the gamma 
radiation observed, the primary CR proton spectrum at the source can be fitted us- 
ing a proton-proton interaction model like the one from Kelner et al. 2006 [ Q) or 
Kamae 2006 lfT3l[l2l . As input for these models, typically a (broken) power-law in 
energy with an exponential cutoff is assumed for the primary CR proton spectrum, 
which forms a spectrum of gamma rays from the decay of neutral pions formed in 
interactions of the CR protons with ambient matter. An example of the shape of the 
primary CR proton spectrum is: 

\£0/ \ ^br/ \ ^cutoff/ 
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Fig. 1 Sketch of the basic idea 



The fit is performed in such a way that the resulting gamma ray spectrum matches 
the observational gamma data. Free parameters of the fit are the spectral indices be- 
low and above the break energy, a/ and a/,, the break energy E^, the cut-off energy 
^cut-off ar, d the normalization of the resulting gamma ray spectrum, ay. A fitted 
gamma ray spectrum for a supernova remnant associated with a molecular cloud is 
shown in figure [2] Since the threshold energy required for the generation of a single 
pion is about 130 MeV/c 2 and a primary CR proton of E p = 1 GeV typically gives 
~ 10% of its energy to a secondary produced, the primary CR proton spectrum ob- 
tained from modelling the GeV-Te V gamma radiation can only be considered known 
down to energies of about 1 GeV. The ionization cross section of molecular hydro- 
gen, on the other hand, has its peak at 100 keV and declines rapidly towards higher 
energies ifTTll . so deriving an ionization rate from the modeled primary CR proton 
spectrum is not trivial. For E p < 1 GeV, the spectrum needs to be extrapolated to 
lower energies carefully, because there is no observational data about this part of the 
spectrum. In a conservative approach, the primary CR proton spectrum is expected 
to change from E~ a to E~ a+i below 1 GeV. Note that the spectral behavior below 
E p < 1 GeV has no significant influence on the resulting gamma ray spectrum due 
to the threshold energy for the production of pions. This modifies the spectral shape 
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Fig. 2 Gamma ray spectrum from pion decay fitted to gamma ray measurements for 
G349.7+0.2. Gamma ray measurements from |8|. Primary CR proton spectral shape 

assumed: jp(E p ) = a p (pj&v) ex p( — rfev)- ^it P arame ters: ^ 2 /dof = 0.88638, 
a y = (820.97822 ± 494.07024)/ (^c) ergMeV- 2 cm" 2 s" 1 , a, = 2.15415 ± 0.11482 . To 
derive « p from ay, see 1181 . 




of the primary CR proton spectrum as: 

te) i(£ ^ £lb) (2) 

(E p >En,), 

where s is varied from +1 .0 to +2.0 , which should offer a lower limit on the spec- 
trum and therefore on the ionization rate. This is due to the fact that an injection 
spectrum of °= p~ a , for a loss term of p « p~ 2 like ionization losses, would 
be modified as °< p 3 ~ a , and the power law indices for the sources discussed are 
s = 1.5 - 2.45 . 

The normalization of the primary CR proton spectrum, a p , is not trivial to calcu- 
late from the resulting gamma ray spectrum normalization ay obtained from fitting. 
In fact, it requires an estimate of the volume of the source, whereas the gamma 
ray spectrum normalization does not depend on this quantity. It enters either in the 
cosmic ray density or in the hydrogen density of the cloud. For a more detailed de- 
scription of the calculation of the primary CR proton spectrum, see Schuppan et al. 

2012 esi. 
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3 Predicting the ionization rate 

Using this spectrum, we can calculate the ionization rate due to primary proton ion- 
ization following Padovani et al. 2009 |17|. Since the spectral behavior of the pri- 
mary CR proton spectrum is not precisely known, the ionization rate is calculated 
for different lower spectral indices s and lower break energies E^. The predicted 
ionization rate only weakly depends on the actual spectral index below the lower 
break, s, but rather strongly on the lower break energy E\\,. This is due to the strong 
decline of the ionization cross section of molecular hydrogen for incident proton 
energies above 100 keV [ 17 1. In our study, we predicted the ionization rates for nine 
molecular clouds known to be associated with supernova remnants. W49B, in par- 
ticular, is expected to show signatures of a largely enhanced ionization rate. For at 
least two objects, the ionization rates are expected to be at least an order of mag- 
nitude above the average of Galactic molecular clouds, aver — 2 x 10~ 16 s _1 
fl6l . In this calculation, there are two quantities left as free parameters. The total 
energy budget in CR protons interacting with the molecular gas, W v , and the hydro- 
gen density of the cloud, «h- The product of these quantities is fixed from fitting 
the resulting gamma ray flux to observational data, W p «h = const., but the val- 
ues of the quantities themselves cannot be derived. Note that this equation can be 
rearranged in such a way that W p «h = PpVdoud n H ~ Pv^cloud — const., where 
p p is the energy density of the cosmic ray protons. Therefore, the calculation can 
be adapted to the quantities for which there are the best estimates available. In our 
calculations, we chose the fixed value of «h = 100 cm~ 3 as standard value for the 
clouds. Should future measurements hint at different densities, the results change as 

This work focuses on the ionization due to primary CR protons, neglecting the 
contribution of primary CR electrons. For a description of the latter, see Padovani 
et al. 2009 ifTTll . Furthermore, the contribution of secondary electrons from primary 
ionization events is neglected here, which typically enhances the primary ionization 
rate by a factor of 1.4 to 2 |9 |. Neglecting these processes, the calculation of the 
ionization rate is performed following Padovani et al. 2009 fTTll as: 

C H2 = / j p (E p )a™(E p )dE p , (3) 

where E m \ B is the minimum energy of protons considered to contribute to the ioniza- 
tion rate and dp " is the ionization cross section of molecular hydrogen for protons. 
£max can be chosen as 1 GeV, increasing it further does not change the resulting 
ionization rate significantly, due to the rapidly declining ionization cross section for 
£p > 100 keV. The value of E m i a , however, is of large importance for the result. 
Generally, protons of kinetic energy E p < 100 keV are not energetic enough to 
penetrate a cloud of «h = 100 cm~ 3 . Due to the rapid decline of the primary CR 
proton spectra used here towards energies E p < £jb, a value of £„,;„ = lOMeVis 
chosen as a conservative approximation [ 10 1 . This lower integration limit becomes 
even more important the lower the value of the lower break energy E\\, is. The latter 
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is varied from 30 MeV to 1 GeV in our calculations. Summed up, the larger the 
lower break energy En, and the spectral index s below this lower break energy, the 
lower the resulting ionization rate. In table [T] the parameters used to calculate the 
ionization rate for all objects considered are shown. The lower break energy, E^,, 
was altered, while the spectral index of the primary CR proton spectrum below the 
lower break energy was chosen as s = +2.0 as the most conservative approach. The 
resulting ionization rate for each object considered depending on the lower break 
energy is shown in figure[3] Changing s from +2.0 to +1.0 does change the results 
by a factor of ~ 1.2 for = 30 MeV up to a factor of ~ 2for£'ib = 1 GeV. This, 
in combination with the value of the lower break energy, is likely to be the largest 
uncertainty in the calculation. As can be seen, for E\\, = 100 MeV, for three objects 
the predicted ionization rate is at least an order of magnitude above Galactic average 
for molecular clouds. This choice of the lower break energy En, might be the most 
reasonable. However, even for the most conservative scenario, i.e. E^ = 1 GeV and 
s = + 2.0, the predicted ionization rate for W49B is orders of magnitude above 
Galactic average for molecular clouds. 



object 


d(k|Xj 


^cloud ( cm ' i ) 


^0 


«/ 


a h ^cutoff 
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Table 1 Table of parameters used to calculate the ionization rate. All values for W p are calculated 
fornn = 100 cm -3 . Gamma ray measurements from: fT1l2l[3ll4ll5l[8l 



4 Ionization signatures 

An enhanced ionization rate of molecular and atomic hydrogen triggers a chemical 
network, forming about 15 to 20 relatively simple molecular ions. In fully molecu- 
lar regions, almost every ionization of molecular hydrogen results in formation of 
an molecular ion. Therefore, the detection of this molecular ion is a very good 
tracer for the ionization rate in a molecular cloud ifTTl , With the ionization rate and 
a few estimates of cloud parameters like the temperature and density, the amounts 
of molecular ions formed in a certain ro-vibrational state can be derived using a 
chemical network and a radiative transfer code as e.g. RADEX lfT9l . In Becker et 
al. 2011 |6|, the first predicted spectra of H^, Ht, and HeH + are shown for a ref- 
erence model cloud. Similar spectra can be predicted for other molecular ions. In 
general, the most reactive ions, such as H^, OH + , and H20 + , will be the most di- 
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Fig. 3 Predicted ionization rate for all objects considered for different lower break energies En, of 
the pirmary CR protons below E p < 1 GeV. 

rect tracers of high ionization. Such a prediction can be made for each molecular 
cloud associated with a supernova remnant. The detection of such relaxation lines 
in spatial correlation with GeV to TeV gamma radiation is a strong hint at hadronic 
origin of both signals. The largest uncertainty in the calculation is the shape of the 
primary CR proton spectrum for energies £ p < 1 GeV. The estimate chosen here 
is rather conservative, so the derived ionization rates are actually lower limits. A 
more precise method to predict ionization rates from CR protons is to propagate the 
CR protons into the cloud, which will be done in future work. Of particular interest 
is the prediction of ionization profiles, since the profiles induced by CR ionization 
are expected to be distinct from those induced by photo-ionization and therefore 
pinpointing the origin of enhanced ionization, at least near the cloud cores, to cos- 
mic rays. With instruments as the Cherenkov Telescope Array and The Atacama 
Large Millimeter/submillimeter Array, the spatial resolution to perform the correla- 
tion study would be available. 
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